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1. Introduction 

Electric dipole moments (EDMs) provide a competitive means to search for new physics 
(NP), complementary to strategies like direct searches at hadron colliders, but also to other in¬ 
direct searches like the flavour-changing processes investigated at the flavour factories. The excep¬ 
tional sensitivity is due to the combination of experimental precision with a tiny Standard Model 
(SM) background. The smallness of the latter is related to the very specific connection between 
flavour and CP violation in the SM,* embodied by the Kobayashi-Maskawa mechanism [1], which 
is very effective in suppressing flavour-changing neutral currents (ECNCs) [2], and even more so 
in suppressing flavour-conserving ones involving CP violation. An exception is provided by the 
gluonic operator 0'qq oc , yielding a potentially very large contribution to hadronic 

EDMs which is, however, strongly bounded experimentally, constituting the strong CP problem. 
To explain its absence, typically symmetries are invoked, involving additional particles. The most 
famous example is the Peccei-Quinn mechanism [3], predicting the presence of axions [4, 5]. 
While these have not yet been observed, it is implicitly assumed in this work when discussing 
hadronic EDMs that the strong CP problem is solved by this or a similar mechanism. The remain¬ 
ing SM contributions then lead to EDMs many orders of magnitude below the present limits, e.g. 
dfj ' ^10 ecTci [6—9]. Importantly, for leptonic EDMs no assumption regarding ^GG is 

necessary; the SM contribution to the electron EDM is estimated to be < lO^^^ecm [10-12]. 
The observation of an EDM with the present experimental precision would therefore clearly con¬ 
stitute a NP signal, especially in the leptonic sector. 

Sakharov’s conditions [13] require the presence of new sources of CP violation to explain the 
observed baryon asymmetry of the universe; while this does not imply sizable EDMs, they are 
generally very sensitive to such sources. In fact, generic NP scenarios usually yield contributions 
that are large compared to experimental limits, implying either a high NP scale or a very specific 
structure for additional CP-violating contributions, similar to the situation in the flavour-changing 
sector. Casting these qualitative statements into reliable bounds on model parameters requires 
knowledge of their relation to the experimental observables - typically (bounds on) frequency 
shifts obtained for composite systems. The calculation of these relations proceeds via a series of 
effective field theories (EETs), see e.g. Refs. [14,15] for recent reviews and references therein. 
Importantly, this approach allows to perform a large part of the analysis model-independently. 
The calculation of the matrix elements of the corresponding effective operators often include large 
uncertainties which have to be taken into account, see Refs. [15,16] for recent detailed discussions. 
Eurthermore, in composite systems different contributions can exhibit cancellations; this issue can 
already be systematically addressed for paramagnetic systems [17,18], and in the future potentially 
also for diamagnetic ones [19]. 

This article proceeds as follows: model-independent constraints are discussed in the following 
section, focusing on paramagnetic systems. In Sec. 3 NP contributions to EDMs are discussed first 
generally, then using Two-Higgs-Doublet models (2HDMs) as a specific example. We conclude in 
Sec. 4. 

*EDMs are T,P-odd; their existence implies also CP violation when assuming CPT to he conserved as we will in 
this article. 
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2. Model-independent constraints from EDM measurements 

The available competitive observables, that is, the EDMs of thorium monoxide (ThO) and 
ytterbium fluoride (YbF) molecules [20,21], thallium (Tl) and mercury (Hg) atoms [22,23] and the 
neutron [24,25] (see also [26]), are related by calculations on the molecular, atomic, nuclear and 
QCD levels to the coefficients of an EFT at a hadronic scale (see, e.g.. Refs. [14,15,27]): 

£%.%., ( 2 . 1 ) 

This operator basis consists of (colour-)EDM operators , the purely gluonic Weinberg operator 
and T- and P-violating four-fermion operators ^ff without derivatives (/^'^ =e,q,q = u,d,s). 
Since these calculations do not depend on the NP model under consideration, this Eagrangian 
is used as the interface between experiment and high-energy calculations: the latter provide the 
model-specific expressions for the Wilson coefficients in Eq. (2.1), with at least one more interme¬ 
diate EFT at the electroweak scale. 

In neutral composite systems, the EDMs of the components are shielded; for non-relativistic, 
point-like constituents this shielding is perfect, therefore measurements for this type of system 
rely on the violation of these assumptions [28]. For paramagnetic systems, relativistic effects 
can actually lead to enhancement factors, if the proton number Z is large enough [29-31], since 
two contributions scale approximately with Z^: these are the ones from the electron EDM and the 
scalar electron-nucleon coupling, Cs-^ Heavy paramagnetic systems can therefore be assumed to be 
completely dominated by these two contributions, allowing a model-independent fit to bound and 
eventually determine both contributions, without the assumption of a vanishing electron-nucleon 
contribution [17,18]. In practice, there are two complications with this approach at present, which 
can however be overcome with additional measurements. Firstly, the ratio of the coefficients of 
the two contributions is necessarily similar for heavy paramagnetic systems [17]. This problem 
can be solved by performing measurements with atoms or molecules with largely different proton 
numbers, such as rubidium and francium atoms. Absent such (competitive) measurements, it is 
possible to assume e.g. the limit from Hg to be saturated by the Cs contribution [18]: this is 
a conservative procedure, since the EDM of this system is typically dominated by colour-EDM 
(cEDM) contributions. While this lead in the past to a similar limit on de like naively setting 
Cs —)■ 0, the unmatched precision of the recent ThO measurement (in the dg-Cs plane) poses the 
second problem, namely that the projections of the 2-dimensional fit remain almost unchanged, 
see Fig. 1 (left). Since this may appear overly conservative, a temporary option is to limit the 
cancellation between the two contributions within the ThO measurement, i.e. to allow the Cs 
contribution to at most saturate the experimental limit n = 1,2,3,... times, see Fig. 1 (right), as a 
compromise between the conservative procedure described previously and the potentially critical 
assumption of Cs = 0. The corresponding fits lead to 

de < 1.0(0.16)10^2’ecm(95%CE), (2.2) 

^Note that Cg depends in general on the considered system. However, for the systems at hand (and more generally 
for heavy paramagnetic systems), it is universal to very good approximation [ 18 ]. 
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using Hg and n = \, respectively. These values should be used when extracting bounds on pa¬ 
rameters from de in any model in which the electron-nucleon contribution cannot be argued to 
be negligible. Furthermore, these 2-dimensional constraints allow to obtain model-independent 



J,/ao~-‘e cm) <4/(10 cm) 


Figure 1; Fit to the recent measurements for paramagnetic systems [20-22], using additionally the Hg 
measurement [23] (blue ellipse on the left) or the assumption of limited cancellations in ThO with n = 1,2,3 
(grey bands, both plots), see text. Plots taken from Ref. [16]. 

constraints on the EDMs of all other heavy paramagnetic systems [18]; examples are 

dRb < 3(0.5) X lO^^^ecm, dc, < 10(2) x < ^0^20) x lO^^e^^j^ (all 95% CL). 

(2.3) 

A significant measurement in one of these systems larger than the first bound would indicate an 
experimental problem, while a value larger than the bound in brackets would indicate large cancel¬ 
lations between the two contributions. These limits are orders of magnitude below existing experi¬ 
mental ones [32-34]. Importantly, present experiments aim at an even better sensitivity [35-39]. 

The extension of this type of fit to all EDM measurements is clearly possible and has been pro¬ 
posed in Ref. [19]. While this is complicated by the many potential contributions - all of the terms 
in Eq. (2.1) are relevant in general -, it should be aimed at in the future. Since model-independent 
bounds/determinations are necessary to determine the specific structure of CP-violating NP contri¬ 
butions - and thereby potentially the model itself -, many measurements in different systems are 
very important. An additional complication for the EDMs of diamagnetic systems and neutrons is 
that the theoretical uncertainties for the relevant matrix elements are often large and can in some 
cases preclude the extraction of conservative limits, for instance on the cEDMs in Eq. (2.1) from 
Hg [16], further motivating complementary measurements. 

3. NP contributions to EDMs 

Reliable limits on parameters in NP models are difficult to achieve. Reasons are, apart from 
the fact discussed previously that presently less competetive measurements than relevant effective 
operators exist, the presence of several contributions to each of these coefficients and the various 
relevant hierarchies, i.e. mass scales, small mixing parameters and loop factors. This complicates 
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semi-model-independent analyses for classes of models and allows strict statements only under 
additional assumptions. 

Generic NP contributions at tree- and one-loop level are in conflict with the stringent exper¬ 
imental limits. On the two-loop level, usually so-called Barr-Zee- and Weinberg diagrams domi¬ 
nate [40^3], which compensate the additional loop factor by avoiding small mass factors. How¬ 
ever, flavour sectors are usually far from “generic”; therefore in some cases also tree-level diagrams 
can be relevant, for example when they involve small mass factors, see below. 

In order to demonstrate these qualitative statements in a specific model, we consider a gen¬ 
eral 2HDM, whose charged current Yukawa couplings can be parametrized in the Higgs basis as 
follows: 


\Pl 

= (3.1) 

where the M, are diagonal mass matrices, W denotes the CKM matrix, and the gp in principle 
arbitrary complex matrices. Below the constraints are given for the elements of these matrices, 
which can be translated into the parameters of any given 2HDM. To be specific and able to relate 
the resulting bounds also to those from other observables, we will furthermore consider the Aligned 
2HDM (A2HDM) [44,45], where the g,- are complex numbers, thereby avoiding FCNCs on tree 
level while still allowing for a rich phenomenology including additional CP-violating phases. The 
couplings of the neutral Higgs mass-eigenstates cp? = {h,H,A} to a fermion species / are denoted 
by yj '; they depend not only again on the Q, but in general additionally on the parameters of the 
Higgs potential, severely complicating the analysis. These couplings fulfill the relations 


£Re 


(p- 

yf 


Im 


y/' 


= ±Im 


(?G/))//(?f (/'))/'/' 


(3.2) 


with a vanishing right-hand side for real g,- (as e.g. the case for 2HDM models) or / = /' 
(more generally fermions of the same family if the g, are family-universal as e.g. in the A2HDM). 
Importantly, the right-hand side is independent of the parameters of the scalar potential. While in 
practical calculations there are mass-dependent weight factors in the sum on the left, the relation 
still holds exactly in two limits: trivially when the neutral scalars are degenerate, but also in the 
decoupling limit [16]. Therefore, in general cancellations can be expected for any mass spectrum 
and the influence of CP violation in the potential is reduced. Clearly, this observation provides 
a protection against large EDMs for models which exhibit new CP-violating parameters in the 
potential, only. Below we will assume this cancellation to occur and evaluate the right-hand side 
with a common weight factor at an intermediate effective neutral Higgs mass M,p. 

In this setup, the situation is typically the one described above: four-quark (tree-level) contri¬ 
butions are subleading, one-loop contributions to (colour-)EDMs are under control (but not nec¬ 
essarily tiny), and two-loop contributions are dominant, but also the tree-level quark-electron cou¬ 
plings are relevant, despite the small mass factors [16,46]. 

The electron EDM receives contributions mostly from Barr-Zee diagrams. The resulting con¬ 
straints on Im(g„ 33 g;*jj) demonstrate the strength of this observable: Im(g* 33 g/n) < 0.05 (for 
n = 1). This already questions the common assumptions that these factors are ^(1). In the A2HDM 
this can be compared to the absolute value of this parameter combination obtained from leptonic 
and semileptonic decays [45,47], which is about a factor 1000 weaker. 
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As mentioned above, also the eonstraint from €$ is relevant: while in this ease the eonstraint 
seems mueh weaker, Im(g^ . 11 ) ~ 15, the bound on the ratio with is again at least a fae- 
tor 100 stronger than the eorresponding one with the eharged-Higgs mass from (semi-)leptonie 
proeesses [45,47]. 

For the neutron, the eonstraint indueed in the eharged-Higgs seetor via the Weinberg operator 
is the dominant one, leading to Im( 5 * 33 ^/ 33 ) < 1. While this does not imply any sizable finetuning, 
it already prohibits large CP-violating effeets in other observables in speeifie models. For instanee, 
while the indireet eonstraint from the branehing ratio mb ^ sy in the A2HDM still allows for a 
sizable CP asymmetry for this proeess, a NP eontribution of |Acp(^ —)• 57 )! <1% follows from the 
EDM bound and the diseussion in Refs. [48,49]. 

These examples show the potential of EDMs, but also their eomplementarity to other searehes, 
sinee only the imaginary parts of parameter eombinations are eonstrained. However, for the eom- 
binations EDMs are sensitive to, they are often the strongen eonstraints available. 

4. Conclusions 

EDMs provide unique eonstraints for the CP-violating seetors of NP models. A potential 
diseovery of a non-vanishing fundamental EDM in any system would be a major aehievement, 
independent of its souree. The interpretation of bounds and potential measurements in terms of 
fundamental theory parameters requires the eareful estimation of theoretieal uneertainties and is 
eomplieated by potential eaneellations on various levels. While this problem ean be addressed 
for heavy paramagnetie systems to extraet the eleetron EDM and sealar eleetron-nueleon eoupling 
model-independently, a similar approaeh ineluding all relevant systems should be aimed at, but 
requires several additional measurements for different systems. 

Eor the oeeuring eombinations of parameters, EDMs typieally provide the most stringent 
eonstraints. We demonstrated this explieitly for general 2HDMs, and more speeifieally for the 
A2HDM, where large CP-violating effeets in other observables are strongly bounded by the exist¬ 
ing EDM limits. Given the present strength of the eonstraints, fortheoming experiments will test a 
erueial part of the parameter spaee and might turn existing bounds into observations. 
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